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We demonstrate that an allometric model for eelgrass leaf-growth rates can be derived from data on leaf 
architecture and growth form. Using this construct, we produced indirect assessments of growth rates of 
leaves that we call projections, which can be easily obtained in terms of allometric parameters and proxy 
values for leaf area, expressed as the product of leaf length and width. These projections of leaf-growth 
rates displayed a high level of correspondence with values observed in our data, as well as with other sets 
of reference data. A comparison with growth rates obtained by using the plastochrone index method 
showed that our model provides more accurate estimations while using a simpler methodology. Our 
results also show that whenever allometric parameters for the scaling of eelgrass leaf dry weight in terms 
of leaf area are available, the proposed model provides an accurate, cost-effective and non-destructive 
alternative to assessments based on traditional or plastochrone methods. 
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1. Introduction 


Constanza et al.. [1] stressed the relevance of ecosystems ser- 
vices as a contribution to human welfare. Based on published data 
for 17 ecosystem services in 16 biomes, they estimated that the va- 
lue of ecosystem services at the level of the whole biosphere. They 
found a lower bound in the range of US$16-54 trillion (10'”) per 
year with an average of US$33 trillion per year. Marine systems 
contributed about 63% of the annual value, with almost half deriv- 
ing from coastal systems. Approximately 25% of the overall contri- 
bution of coastal ecosystems corresponded to algal beds and 
seagrasses, which are submerged angiosperms that are considered 
a characteristic feature of shallow environments in tropical and 
temperate zones. 

Petersen [2,3] pioneered research aimed at characterizing the 
ecological relevance of seagrass species. His studies provided 
groundbreaking perspectives on the functional attributes of seag- 
rasses, emphasizing their role as vigorous primary producers. 
Nonetheless, half a century passed before marine scientists fully 
acknowledged that seagrass meadows are among the most produc- 
tive ecosystems, rivaling tropical agriculture [4-10]. Although high 
primary productivity is a key feature in the trophic role of seagrass 
meadows [11-13], they also furnish an essential structural ele- 
ment by providing habitat for the growth of epiphytes [14,15] 
and shelter for animal communities [16,17]. Seagrasses are thus 
not only useful indicators of the health status of coastal zones, 
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but also their ability to modify their own environment, particularly 
in relation to the processing of sediments, allows them to be used 
in coastal management strategies aimed at preserving or improv- 
ing environmental quality. It is this important attribute that has 
granted them the title of ecosystem engineers [18-21]. Neverthe- 
less, it is their contribution to nutrient cycling within estuaries 
which represents the largest proportion of the economic benefit 
performed by seagrass meadows [22-24,1]. 

Eelgrass (Zostera marina L.) is a particularly appropriate species 
for the study of seagrass ecology, not only because it is the domi- 
nant species along the coasts of both the North Pacific and North 
Atlantic [25], but also because eelgrass communities are recog- 
nized as among the richest and most varied in the abundance of 
sea life [26]. Indeed, this cosmopolitan macrophyte was found to 
produce up to 64% of the total primary production of an estuarine 
system [27]. Therefore, the assessment of eelgrass biomass vari- 
ability through time provides valuable information on the dynam- 
ics of both its trophic and structural roles. The basic unit for 
studying biomass production in eelgrass is the shoot (see Fig. 1), 
which includes leaves, rhizomes and roots. A rhizome is a horizon- 
tal stem that grows along the sediment, forming roots as it elon- 
gates. Root emergence occurs at leaf scars, also known as 
rhizome nodes. Leaves appear from a meristem lying within a pro- 
tective sheath. This structure also holds the leaves together and is 
attached at the actively growing end of the rhizome. The produc- 
tion of leaves and rhizome nodes is linked in such a way that every 
leaf produced is associated with a rhizome node. Hence, the overall 
production of shoots can be estimated by measuring the produc- 
tion of leaves [28], which makes the determination of the growth 
rate of leaves a fundamental entry for the assessment of eelgrass 
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Fig. 1. A pictorial representation of an eelgrass shoot. The segmented stem 
portrayed is called a rhizome. Each segment is bounded by nodes also called 
rhizome scars. These nodes are produced when new leaves are formed. The bundle 
of leaves is attached to the sheath, a specialized structure located at the end of the 
rhizome which contains a meristem where new leaves are produced. This drawing 
was adapted from [63]. 


populations. The rate of production of leaf biomass in eelgrass is an 
indicator variable for environmental influences on growth and also 
on the overall functions and values of this important seagrass spe- 
cies [29]. 

The forcing by environmental variables on the dynamics of Z. 
marina has been well documented in the literature [30-32]. Light 
availability, temperature and dissolved nutrients have been 
claimed as the most important variables for explaining the ob- 
served variability [29,33]. In particular, sea-surface temperature 
drives changes in aboveground biomass, and is highly correlated 
with leaf dynamics [29,34-36]. Hence even when light and nutri- 
ents are not limiting, temperatures beyond the upper limit of its 
tolerance range can provoke severe negative effects on eelgrass 
growth [37]. Global warming is expected to result in warming 
and rising seas, thus reducing the availability of both light and 
nutrients underwater [34]. Indeed, the onset of warm ENSO events 
has been shown to dramatically diminish eelgrass growth [37]. An- 
other concern for the health of eelgrass populations is associated 
with deleterious anthropogenic influences. The loss of habitat for 
this macrophyte has been noted worldwide, with catastrophic 
losses within the past few decades [38-40]. These concerns have 
driven efforts aimed at preserving their populations, and several 
workers have relied on transplantation projects as a way to restore 
lost habitats [41-43,25]. 

Due to seasonality effects, the study of leaf growth in Z. marina 
demands extensive sampling during the entire yearly cycle. Tradi- 


tional assessment methods require destructive and time consum- 
ing dry-weight determinations. The destruction of shoots during 
early stages of growth might provoke undesirable alterations in 
transplanted populations. To eliminate these inconveniences, we 
present a non-destructive and cost-effective method for the accu- 
rate estimation of growth rates, which uses an allometric model 
for the representation of leaf biomass in terms of the product of 
leaf length and width. The use of this construct could make the re- 
moval of shoots and the traditional dry-weight determinations 
unnecessary, and thus greatly simplifying conventional leaf- 
growth assessments including the plastochrone method [30,44]. 


2. Data and related variables 


The primary data set used in this study was collected form April 
1998 to April 1999 in a Z. marina meadow in the Punta Banda estu- 
ary, located near Ensenada, Baja California, Mexico. For a detailed 
description see [35]. Using the technique of Kentula and McIntire 
[45], we marked approximately 40 shoots at each sampling time 
t. At a date t+ At, with At fixed at 2 weeks, we removed the 
marked shoots that were still present and marked a new sample. 
In what follows we will let I(t), h(t), a(t) and w(t), respectively, 
denote length, maximum width, area and dry weight of a typical 
Z. marina leaf. For each one of the leaves on a recovered shoot, 
the corresponding I(t + At) values were directly obtained by mea- 
suring the distance between the ligula at the top of the sheath and 
the leaf tip. We also measured the leaf elongation increments 
Al = I(t + At) — I(t) which are given by the displacement of the 
marking point from the ligula in the technique we used. Observa- 
tions indicate that a typical Z. marina leaf has a long and narrow 
belted shape. This architecture can be explained by assuming that 
leaves grow initially both in a longitudinal and a transverse dimen- 
sion, but while width quickly reaches a maximum value, length 
continues to increase towards an asymptotic value [29]. For a given 
maximum width, a great variability in corresponding leaf lengths 
has been observed. Hence we judged that measurements taken at 
a point halfway between the leaf tip and the ligula provided repre- 
sentative estimations for h(t + At). This is consistent with other 
estimations of eelgrass leaf width reported in the literature, e.g., 
[46-48]. Following a customary approach, e.g., [49-54] we choose 
the product of I(t + At) and h(t + At) as a reasonable proxy value 
a,(t + At) for the real leaf-area value a(t + At). Moreover, for the 
biweekly marking interval considered, we also observed that an 
assumption of negligible values for Ah is fairly reasonable. Hence, 
we regarded the product h(t + At)Al as a proxy value Aa, for the 
real leaf-area increment Aa = a(t + At) — a(t). The leaf biomass 
increments Aw=w(t+At)—w(t) associated with a leaf-area 
increment Aa were estimated through conventional procedures 
[44]. The observed values for Aw on a given shoot s were used to 
estimate the related mean rate of leaf growth over the interval 
(t,t + At) These quantities are denoted here by means of the sym- 
bol L..(t, At) and formally defined through, 


;AW 


Isg(t, At) = =A, (1) 





where 5°, indicates summation over a given shoot. This notation 
convention, along with the use of a s subscript to label a generic 
shoot will be kept throughout. 

The number of leaves produced by a shoot after it was marked 
was used to calculate the plastochrone interval using the formula 
in Jacobs [30]. These values were used to obtain alternative estima- 
tions for L,,.(t, At) according to the plastochrone method described 
in Gaeckle and Short [44]. Finally, for our own and the plastoch- 
rone method determinations, we calculated the corresponding 
overall rates of leaf growth. These are referred to here by means 
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of the symbol L,(t, At), and obtained averaging the L,.(t, At) values 
given by Eq. (1) over the number of collected shoots, i.e., 


L,(t, At) = | , 


where NS(t, At) stands for the number of collected shoots. 

In addition to the data presented in this study, we also used 
three other independent data sets for comparison purposes. The 
first, reported in Solana-Arellano et al. [55], was collected in San 
Quintin Bay during 1992-1993 and included lengths, widths and 
dry weights of complete leaves. A second data set, reported in 
Solana-Arellano [56] presented leaf lengths, widths and dry 
weights that were also obtained from our Punta Banda study site 
in 1999. Finally a third data set collected in 1987 in San Quintin 
Bay was obtained from Ibarra-Obando [57]. This includes plastoch- 
rone interval values, leaf lengths, widths, estimations for Aa and Al, 
as well as observed values for corresponding L,(t, At) rates. 


3. Theoretical methods 


Following Solana-Arellano et al. [58,59] we can show (see Eq. 
(A2)) that w(t) can be scaled in terms of a,(t) through a model of 
the form, 


w(t) = pap(t)*h(t)’. (3) 


where p > 0, and « and @ are constants. It can also be shown (see 
Eqs. (A2)-(A5)) that, on the basis of the observed leaf architecture 
and growth form, we can derive from the above equation an allome- 
tric model for w(t) expressed solely in terms of a(t). This is, 


w(t) = Cap(t)", (4) 


where c > 0 and m are constants. Using this result, we can obtain 
indirect estimations for the corresponding leaf dry-weight incre- 
ments Aw produced over the interval (t,t + At). In fact, in a straight- 
forward fashion from Eq. (4) we can obtain for Aw 


Aw = Cap(t + At)'"5p(t, At), (5) 


where the weighting factor 6,(t, At) is given by 


5, (t, At) (1 (1 ates) ): (6) 


The a,(t + At) values can be directly obtained using measure- 
ments of leaf length and width on shoots collected. Since direct 
estimations for a,(t) are unavailable, we must use leaf-marking 
data to produce these values. Indeed, the measured proxy values 
Aa, will determine a,(t) through the _ relationship 
y(t) = dp(t + At) — Aap. 

Adding the projected increments Aw on a given shoot and divid- 
ing by At we obtain projected values L,;.(t, At) for the observed 
mean leaf-growth rates L.,.(t, At) defined by Eq. (1) formally, 


di sCAp (t a3 oe At) . (7) 


Averaging these values over the number of collected shoots we 
obtain a projected value Ly,(t, At) for the overall mean leaf growth 
L,(t, At) defined by Eq. (2), e.g., 








Lysg(t, At) 








It is worth pointing out that the reliability of the estimation gi- 
ven above depends on the consistency of the fitting of Eq. (4) to 
eelgrass data. 

In Appendix we also show that there are suitable allometric 
constants b and k such that the result of Eq. (4) can be also 


expressed directly in terms of the real leaf-area values a(t), i.e., 
we also have 


w = batt)", (9) 


for b > 0 and k constants. Moreover, using the symbol Lgse(t, At) to 
denote projected leaf-growth rates expressed in terms of a(t), we 
obtain (see Eqs. (A6)-(A12)), 


So, ba(t + At)*5q(t, At) 








Losg(t, At) = , (10) 
where 

: 7 a(t) 
ba(t, At) = ( i). (11) 


Moreover, this last estimation device does not require the use of 
leaf-marking techniques, which means a further simplification to 
the projection method given by Eq. (7). This last estimation device 
can be used whenever an automatic leaf-area meter is available. In 
such case we can produce direct and highly accurate estimations 
for a(t) [49]. 


4. Results 


We fitted Eq. (3) to the Solana-Arellano et al.[55] data to obtain 
p = 0.000006, «= 1.24, 0 = —0.44 with a determination coeffi- 
cients of R* = 0.93. The fit for the present data set produced 
p = 0.000004, « = 1.36 and @ = —0.61 and R* = 0.87. Meanwhile, 
for the data reported in Solana-Arellano [56] we obtained 
p = 0.000004 and 0 = -0.57, «= 1.44 with R? = 0.95. Residual 
analyses for all these fittings showed that the residuals behave in 
accordance with the assumption of a general regression model, 
thus indicating that our model is consistent. 

Using the above data sets, we also corroborated that eelgrass leaf 
architecture upholds a consistent identification of the simpler 
allometric model given by Eq. (4). The values for the parameters 
involved were c= 0.000004 (std-err.=0.02) and m=1.21 
(std-err. = 0.003) with R? = 0.95 for the Solana-Arellano et al. [55] 
records. For the present data set we found c=0.000003 
(std-err. = 0.002) and m = 1.25 (std-err. = 0.16) with a determina- 
tion coefficients of R? = 0.92. Meanwhile for the Solana-Arellano 
[56] data we found c=0.0000025 (std-err. = 0.019), 
m = 1.3 (std-err. = 0.15) and R? = 0.96. For all these fittings, the 
analysis of residuals shows consistency. The present data set and 
the Solana-Arellano et al. [55] are from different sites, so we judged 
it more informative to include in (Fig. 2) only the corresponding 
fittings. 

Because the data for this study and that of Solana-Arellano [56] 
were obtained at the same study sites, we tested for differences in 
the allometric parameters c and m through time, and found no sig- 
nificant differences (p = 0.49 for c and p = 0.71 for m). This sug- 
gests that, for our study site, a single representative sample 
including w(t) and a,(t) values, will produce estimations of the 
parameters c and m that are expected to be valid through a com- 
plete yearly cycle of eelgrass growth. Moreover this property of 
time invariance in c and m might also be extended to data taken 
from a different study site. In fact, a test for statistical differences 
of fitted c and m values between the Solana-Arellano et al. [55] 
and the Solana-Arellano [56] data set also resulted in non-signifi- 
cant differences in allometric parameters between San Quintin 
Bay and the Punta Banda estuary (p = 0.50 for c and p = 0.62 for 
m). The time span between these studies was 7 years. During these 
studies, different environmental conditions prevailed at the host 
sites (http://www.cdc.noaa.gov.people.klauss.wolter/MEI; http// 
ferret.pmel.noaa.gov/NVODS/servlets/datasets), but only negligible 
differences were noted in the fitted parameters. This indicates that 
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Fig. 2. Fittings of the allometric model w(t) = ca,(t)” (a) for the Solana-Arellano et al. [55] data (c = 0.000004, m = 1.21 and R’ = 0.95) and (b) for the present study data 


(c = 0.000003, m = 1.25 and R? = 92). 


the assumption of time invariance for c and m in our study site 
could be extended to San Quintin Bay (which is approximately 
100 km to the south). Moreover, we can also assume that parame- 
ters fitted to data sets taken at either of these sites - even during 
different years — could be used to project leaf-growth rates, for 
example in cases where the data required for the fitting of the 
model of Eq. (4) are missing. This criterion was used to produce 
projected values for leaf-growth rate from the data reported in 
Ibarra-Obando [57]. 

Once the parameters c and m were identified, we proceeded to 
apply our projection method (cf. Eq. (8)). We recall that for the 
Solana-Arellano [56] and Solana-Arellano et al. [55] data sets we 
could not obtain projected values Lp.(t, At) due to lack of estimated 
values for Aw and Aa. Fig. 3 displays a comparison of projected 
Lyg(t, At) and observed leaf-growth rates L,(t) for the present and 
the Ibarra-Obando [57] data sets. For the latter data set we relied 
on parameter values c and m fitted from Solana-Arellano et al. 
[55]. The level of correspondence between projected and observed 


values shown in Fig. 3(a) and (b) indicates that the results were 
trustworthy. This level of correspondence is remarkable if we com- 
pare with alternative assessments obtained by means of the named 
plastochrone method [30,44] (see Fig. 3(c) and (d)). 

The robustness of our method was also explored through cross 
validation; that is, by interchanging fitted parameters. We pro- 
jected the observed leaf-growth rates for the Punta Banda estuary 
using the parameters c and m fitted from the Solana-Arellano [56] 
data and compared those with the observed values reported here 
(Fig. 4(a)). Similarly we compared the above referred projections 
with those produced by our own data set (Fig. 4(b)). In the same 
manner, using c and m found for the San Quintin Bay data and re- 
ported in Solana-Arellano et al. [55] we projected leaf-growth rates 
for Punta Banda estuary and compared with projections obtained 
using the present data (Fig. 4(c)). Finally we projected leaf-growth 
rates for San Quintin Bay using c and m obtained from Punta Banda 
estuary data and compared with those projected using c and m fit- 
ted for the Solana Arellano et al. [55] data set (Fig. 4(d)). For an 
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Fig. 3. Plots of projected (dashed line) and observed (continuous line) mean leaf-growth rates through time. The plots on (a) and (b) correspond to values projected using Eq. 
(8). Panel (a) corresponds to the present data set. Panel (b) shows projections obtained using the Ibarra-Obando [57] data set. The plots on (c) and (d) represent mean leaf- 
growth rates obtained by means of the plastochrone method described in Jacobs [30] and Gaeckel and Short [44]. Panel (c) shows estimations for the present data set. Panel 


(d) displays assessments corresponding to the Ibarra-Obando [57] data. 


additional strength criterion, we tested for differences between 
observed and projected mean annual growth rates at each study 
site and for all projections, and found no significant differences 
(p > 0.09 for all). This seems to indicate that, at each site, the ob- 
served mean annual growth rate is consistently represented by 
the corresponding projected values. 


5. Discussion 


The advantages of allometric relationships in seagrass research 
have been pointed out by several authors, e.g., [30,6,61,46,35]. In 
particular, for several pairs of representative X and Y seagrass vari- 
ables, Duarte [62] found consistent representations of the form 
Y = «X? with a and £ constants. This study shows that eelgrass leaf 
architecture sustains the identification of a similar scaling relation- 
ship w(t) = ba(t)* for leaf dry weight w(t) in terms of the area a(t). 
Moreover, leaf-growth rates Lp,(t,At) derived from the model 
w(t) = ca,(t)” displayed significant levels of correspondence with 
observed values. When compared with pertinent assessments ob- 
tained by means of the plastochrone method, the proposed allome- 
tric approach turned out to be of greater accuracy. And the result of 
a test of non-significant differences for c and m using the refer- 
enced data sets suggests that these allometric parameters can be 
considered as time invariant. This means that a single representa- 
tive sample of leaf areas and corresponding dry weights can be 
used to produce projections which are be expected to be valid 
through a complete observation period. Moreover, we corrobo- 


rated that parameter values can be interchanged among the sites 
considered, while still obtaining robust estimations. This suggests 
that the proposed allometric method can provide a simpler and 
more reliable non-destructive alternative for assessments of eel- 
grass leaf-growth rates. Because the estimation of leaf biomass 
production rates is an important input for assessing the status of 
eelgrass populations, the proposed method represents an added 
benefit of scaling relationships in eelgrass research. For example, 
in transplanted areas, determinations of leaf-growth rate are a fun- 
damental entry to assess the restoration of ecological functions. A 
complete description of growth dynamics must consider estima- 
tions taken over the early stages of the transplant experiment, 
making it essential to rely on non-destructive methods. It is in this 
regard that our projection method could be most useful. In any 
case, our approach amounts to a reliable and cost-effective alterna- 
tive to traditional or plastochrone-based determinations. More- 
over, the L,,(t, At) projections for the L,(t, At) rates can be easily 
obtained, because the parameters c and m are derived from stan- 
dard regression procedures and the a,(t) values calculated from di- 
rect measurements of leaf lengths and widths and leaf-marking 
data. 

The result of Eq. (9) shows that we can use other estimations of 
leaf area, different from a,(t), such as those produced by using an 
automatic leaf-area meter [48]. These estimations for a(t) are faster 
to obtain and may produce smaller approximation errors than 
those associated with a,(t). This procedure could also improve 
the quality of the fit of Eq. (9) over the one in Eq. (4), thereby 
resulting in more accurate projections. Because periods of low-tide 
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Fig. 4. Comparison of leaf-growth rates projected by means of Eq. (8) interchanging the parameters c and m. (a) Observed rates (dashed lines) and projected ones for the 
Punta Banda estuary using the Solana Arellano [56] data (continuous lines). (b) Rates projected for the Solana Arellano [56] data (dashed lines) and projections using the data 
set for the present study (continuous lines). (c) Rates projected with the present study data (dashed lines) and those projected for the same site but using the parameters c and 
m fitted for San Quintin Bay (i.e., using the Solana Arellano et al. [55] data (continuous lines). (d) Rates projected for San Quintin Bay using data in Ibarra-Obando [57] with the 
parameters c and m fitted for Punta Banda (dashed lines) and projections for San Quintin using the Ibarra-Obando [57] data with the parameters c and m found from the 


fitting of the Solana Arellano et al. [55] data set (continuous line). 


exposure are short, and plants may be submerged, it might be 
difficult to obtain a suitable number of replicates for a(t) through 
in situ measurements. Under these circumstances, one could rely 
solely on the leaf-marking approach described here. Even in the 
case in which one is limited to the use of leaf marking, our results 
show that reliable projected values for the L(t, At) rates can be ob- 
tained without the tedious dry-weight determinations of tradi- 
tional or plastochrone assessment methods. The proposed model, 
moreover, can be applied to other seagrass species that have leaf 
architecture similar to eelgrass in a straightforward manner. 

The results presented here raise the possibility that allometric 
parameters obtained at different sites and under different envi- 
ronmental conditions could be used in our paradigm to create 
independent and reliable projections of leaf-growth rates, thus 
amounting to a general, robust and non-destructive assessment 
method. The fact that the referenced data sets produce similar 
estimates for c and m suggest that this could be possible. Never- 
theless, a generalization from the present results amounts to 
considering that c and m are species-specific properties and this 
must necessarily be substantiated by a larger data base, covering 
a larger range of distribution of eelgrass. This will be the subject 
of further research. 
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Appendix 


Solana-Arellano et al. [58,59] generalized the model of Ham- 
burg and Homann [46] and used independent data sets to corrob- 
orate that w(t) can be scaled by using an allometric model of the 
form, 


w(t) = pl(t)"h(e)’, 


where p > 0, « and f are constants. We will show through that we 
can also express the above bivariate scaling of w(t) solely in terms 
of the considered proxy area a,(t). In order to obtain such a repre- 
sentation we first notice that Eq. (A1) takes the equivalent form, 


w(t) = pap(t)*A(t)’, 


where 0 = f — a. Now, let’s consider a particular marked leaf. Then 
we will have w(t) > 0, and direct algebraic manipulation of Eq. (A2) 
yields, 


st (ot) (ay 


Since we observed that during the chosen marking interval eel- 
grass leaf architecture and growth form resulting in negligible val- 
ues for Ah, we might expect the ratio involving a,(t) in Eq. (A3) to 
be dominant. On the other hand, if we assume that w(t) can be 
allometrically scaled in terms of a,(t), then, there will be a constant 
m [60] such that, 


(Al) 


(A2) 








(A3) 
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1dw_ mda, 
wdt ay dt’ 
and from this equation we obtain, 


w(t + At)  dp(t + At)” 
w(t) a,(t)” 








(A4) 


Hence a dominance of the ratio involving ap(t) in Eq. (A3) is con- 
sistent with the assumption of an allometric scaling of w(t) in 
terms of a,(t). Furthermore, assuming that Eq. (A4) holds for all t 
values, then both sides must be equal to a constant. Let c be such 
a constant. Then we will have 


w(t) = ca,(t). (A5) 


The result for Aw in the form of Eq. (A5) can be also obtained for 
leaf-area values a(t). In fact, recalling that a,(t) is an estimator for 
a(t) there must exist an approximation error &,(t) such that, 


a(t) = a,(t) + €)(t). (A6 
Moreover for the ratio 
s(t) = ca (A7 


an additional consistency condition 


s(t) <1, (A8 





is required. 
Solving for a,(t) in Eq. (AG) and replacing in Eq. (A5), we can 
obtain, 


w(t+ At) /a(t+ At)\™ s(t + At) —s(t)\” 
Tce aCe A me cme ve 








By virtue of the consistency condition in (A8), the right-hand 
side of Eq. (A9) will be dominated by the ratio a(t + At)"a(t)"”. 
Hence an elaboration such as the one used to establish Eq. (A5) will 
pertain, so we can also propose a general model, for eelgrass leaf 
dry weight w(t) and corresponding leaf-area values a(t), e.g., 


w(t) = ba(t)*, (A10) 


where b and k are constants. In a manner like that established for 
Eq. (7) we can obtain an equation for the representation of pro- 
jected rates Las, (t, At) in terms of a(t) and the allometric parameters 
b and k, this will be, 


<ba(t + At)*5a(t, At) 





Losg(t, At) = » 


x (A11) 
where 
k 
da(t, At) = (1- =a (A12) 
a(t + At) 
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